DYNAMICAL ASPECTS OF MEAN FIELD PLANE ROTATORS 
AND THE KURAMOTO MODEL 



LORENZO BERTINI, GIAMBATTISTA GIACOMIN, AND KHASHAYAR PAKDAMAN 

Abstract. The Kuramoto model has been introduced in order to describe synchroniza- 
tion phenomena observed in groups of cells, individuals, circuits, etc... We look at the 
Kuramoto model with white noise forces: in mathematical terms it is a set of N oscilla- 
tors, each driven by an independent Brownian motion with a constant drift, that is each 
oscillator has its own frequency, which, in general, changes from one oscillator to another 
(these frequencies are usually taken to be random and they may be viewed as a quenched 
disorder). The interactions between oscillators are of long range type (mean field). We 
review some results on the Kuramoto model from a statistical mechanics standpoint: we 
give in particular necessary and sufficient conditions for reversibility and we point out 
a formal analogy, in the N — > oo limit, with local mean field models with conservative 
dynamics (an analogy that is exploited to identify in particular a Lyapunov functional in 
the reversible set-up). We then focus on the reversible Kuramoto model with sinusoidal 
interactions in the N — * oo limit and analyze the stability of the non-trivial stationary 
profiles arising when the interaction parameter K is larger than its critical value K c . We 
provide an analysis of the linear operator describing the time evolution in a neighbor- 
hood of the synchronized profile: we exhibit a Hilbert space in which this operator has a 
self-adjoint extension and we establish, as our main result, a spectral gap inequality for 
every K > K c . 
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1. Introduction 

The Kuramoto model {e.g. [HET]) is denned by the set of coupled stochastic differential 
equations: 

K N 

d</#) = & dt - - J2 sin ~ Pi (*)) dt + ° dw j(t) > (1-1) 

i=l 

for j = 1, 2, . . . ,N, where 

(1) {wj(-)}j = i } ... s N is a family of independent and identically distributed standard 
Brownian motions. We refer to this source of randomness as thermal noise. 

(2) £ = {S,j}j=i,...,N is a family of independent identically distributed random variables. 
This is another source of noise, and we refer to it as disorder. 

(3) K is a real parameter and a > 0. 

We stress from now that we consider the stochastic evolution (jl.ip once a realization 
of the disorder variables £ is chosen, so the disorder is of quenched type. Moreover the 
law of {^ (•)}j=i,...,iV and of the initial condition (specified below) does not depend on the 
values of the disorder variables. 
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Remark 1.1. It will be at times interesting to discuss the role of the sine drift in the 
model. We will therefore refer to a /i-model when iCsin(-) is replaced by a smooth, i.e. 
C°°, 2-7r-periodic function h{-). 

The variables are actually angles, so we focus on 9j|(i)mod(2-7r), which is an element 
of § := R/(27rZ). The existence and uniqueness of a unique (strong) solution to the system 
(jl.ip . when the initial condition {v|(0)}j=i v .. 5 Ar G M. N and {w^-)}.,-^...^ are independent 
and {Vj(0)}j=i,...,Ar are square integrable random variables, is a standard result. In our 
case we may therefore choose {<p^j(0)}j=i,...,N arbitrarily distributed provided that it is 
concentrated on [0,2-71")^. 

The main result of this work is on the model in which there is no disorder, that is the 
case in which the law of £i is degenerate, so that = £ for every j, with £ is a real 
constant. In this case, with the change of variable <Pj(t) := ipj(t) — £i we have 

K N 

= --5>n((^)-^(t)) dt + ad Wj (t). (1.2) 
i=l 

Disregarding the disorder is actually a major simplification first of all because, if a > 0, 
the system (jl.2p is reversible with respect to the (Gibbs) probability measure 

fJ,N,K(d<p) := — — exp ( -^-H N (ip) ) Ajv(d^) , (1.3) 



where tp G S , A^v is the uniform probability measure on (that is the iV-fold product 
of Lebesgue measures normalized by (2^)^), 

^ N N 

H N ((p) := ~^J2J2 cos ^ ~ fi) ' ( L4 ) 

i=l j=l 

and Zjv,k := J s n ex.p(—2Ka~ 2 H at (</?)) Aat( dip) is the partition function. In fact, the 
generator Lk,n of the dynamics (|1.2|) acts on twice differentiable functions i 7 : S N — > R 
as 

L K}N F(ip) = — y—^~2 — a a ' ( L5 ) 

2 f=! d $ i=l 9 ^ 

and one directly verifies the symmetry J F Lk,nG dfix,N = J G Lk,nF d(J-K,N for F,G £ 
C 2 , which implies that /J-k N is invariant for the dynamics [20^ I13j. 

The measure Hk,n is the Gibbs measure of a classical statistical mechanics model: the 
mean field spin XY model with single spin state space S, i.e. mean field plane rotators 

EES. 



Remark 1.2. It is important to notice that the system (jl.ip is not reversible unless 
£i = for every i. Even the case = const. ^ for every i is not reversible, but, as 
we have argued, it maps to a reversible system. Notice in fact that, unless £ = 0, the 
transformation ip^ i— ► cpj — £i maps to a system with time dependent interactions. This 
strongly hints to the absence of reversibility and it is indeed the case, but proving such a 
statement is more delicate: we address this point in Section [J] below. The aspect that we 
want to stress here is the disorder induced non-equilibrium character of the full Kuramoto 
model. 
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1.1. Empirical measure and the large N limit. Since we focus on the a > case, 
there is no loss in generality in choosing a = 1 and we will do so from now on. We 
introduce the empirical measure 



v N,t\ 



1 - 



(1.6) 



3=1 



and observe that, by Ito's rule, for every F S C 2 (§) and t > 



s 



F(e)u N>t (de) - / F(d)v N>0 (M) 



-K f f F'(9) S m(6-6')v N , s {d6)v N , s (d6')ds + l [ [ F"(9)u N>s ( d9) ds + M N>F (t) , 
Jo Js 2 * Jo Js 

(1.7) 

where Mn,f(') is a continuous martingale with quadratic variation at time t, (Mjv,f)(^), 
equal to iV _1 f Q f g (F'(9)) 2 i>N >s ( d0) ds. Therefore, by Doob's inequality, for every T > 
we have that E[sup te[0jT ] (M N)F (t)) 2 ] is bounded by {M N)F )(T) < T\\F' / N : this guar- 
antees that the thermal noise disappears as ./V — > oo, so that the limit of the empiri- 
cal measure, if it exists, is not random. To make this precise we introduce the space 
C°([0,T]; 7V4i(§)), where A^i(S) are the probability measures on § equipped with the 
topology of the weak convergence, and observe that if a subsequence of {i^-jjveN (of el- 
ements of C°([0, T]; A^i(S))) converges to a limit v., we have that for t 6 (0, T] and every 
F G C 2 (§) 



F(0)vt(M)- / F(^ (d#) = 

\! f F"(9)u s (d9)ds -K I I F'{0)sm{9-9')v s {d9)is s (d9')ds. (U 
2 Jo js Jo Js 2 

This is a weak form of the equation 



dM6) = 2^9^ + K d9 



S m(9-9')q t {9')d9')q t {9) 



(1.9) 



when fj(d#) = <ft(#) d#. So that if we assume that vn,o converges to a non random limit 
and if there is a unique solution to (jl.8p , then the evolution is non random and determined 
by (THJ. 

More precisely, we have the following: 

Proposition 1.3. If there exists uq £ A^i(S) such that for every e > and every F £ 
C°(§;R) we have 



lim ] 



F(9)v Nfi (dO)- I F(0)u o (d9) 

s 



then for every t > we have that for every £ and F 



lim I 



F{9)u N)t {d9) - I F(9)u t (d9) 



> £ 



> £ 



0. 



0. 



;i.io) 



(i.ii) 



where v. is the unique solution of (|1.8|) . Moreover, for every t > the measure v% is 
absolutely continuous with respect to the Lebesgue measure with (strictly) positive density 
qt(-) and the function (i, 9) qt{9), from (0, oo) x § to (0, oo), is smooth and solves 
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Proposition 11.31 is a particular (and particularly easy) case of far more general results 
(see for example [8j US])- The derivation goes along proving tightness of {^a^Ia^n and 
then proving uniqueness for the limiting equation (|1.8|) . In our case such an equation is 
particularly nice and the evolution is smoothing, so that even if the initial datum is not a 
function (i.e vq is not absolutely continuous with respect to the Lebesgue measure) or it is 
not smooth, vt(d9) = qt{9) d8 and qt(-) G C°°(S) for every t > 0. These analytic aspects 
are taken up with more details in Section [3] (Proposition 13, 1|) 



Remark 1.4. It is however important to recall here that Proposition ll.3l can be generalized 
to cover the disordered case (|1.1|) . We refer to [3 [10] for precise statements, but, roughly, 
if the law of the random variable £i is denoted by /x(d£) (let us for example assume that 
£i is bounded), the empirical average at time t > converges as N — * oo to a measure 
with density j R qt(9; £)/z( d£), where {qt(9; £)}t>o,0eS,£eR is the unique solution to 

a fa t\ l 9 ^t(0\0 , d 

« m = ^ 



jf ^Ksm(9-e , )q t (0';e)de'\ /i(d£') + ^ <ft(0;£) 



d0 ' 

(1.12) 

for every £ in the support of /i. We have of course assumed, for simplicity, that uq is 
absolutely continuous with respect to the Lebesgue measure. 



Remark 1.5. The Kuramoto limit evolution (jl.9p comes up also as mesoscopic scaling 
limit for the density of Kac models with conservative dynamics [9] , when the interaction 
potential is chosen equal to cos(-). This is just a formal analogy, but it is going to be crucial 
for the sequel. 

1.2. Stationary profiles. By the regularizing character of the evolution (Proposition l3.ip . 
the stationary solutions to (|1.8j) coincide with the stationary solutions to (|1.9h (we are of 
course interested only in non-negative solutions of total mass equal to one: Proposition l3.ll 
guarantees also the positivity of stationary solutions). Let us notice moreover that if q(-) 
is a stationary solution, then q{- + 6q) is a stationary solution too, for arbitrary choice 
of 6q. This is due to the invariance of fjl .21) under rotations (that can of course be read 
also out of (jl.4p ). Note that q(-) = l/(2n) is a solution to (II. 2h . regardless of the value 
of K but there may be more solutions: in fact every stationary solution can be written as 
<?(• + $0 ) f° r some 9q E [0, 2n) and 

•Ym = exp(21fr cos(fl) 

q{ ' ■ f s exp(2Krcos(e')dO' , [ ' 

with r a non-negative solution to 

, , L cos(0) exp(x cos(0)) d6 
r:=V(2Kr), with V(x) := Js K > PK | » . 1.14 

J s exp(x cos(c/JJ dti 

In general, there is more than one solution to f)l . 14|) : in fact, there can be at most two, 
more precisely there is only the trivial solution r = for K < 1 and there is also a second 
solution r > if K > 1. This is because ^'(0) = 1 and because : [0, oo) — > [0,1) 
is strictly concave [16]. In terms of stationary solutions, this means that for K < 1 only 
the flat (incoherent) profile l/(27r) is stationary, while for K > 1 also {g(- + 0o)}6» o eS 
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is a family of stationary solutions (they are the solutions that exhibit the coherence or 
synchronization of the system). 

The result we just stated, that is (|1.13|) - (|1.14|) . is a classical one in the sense that it is 
of course closely linked to the solution of the mean field planar rotator model [18] (result 
completed by the concavity result proven in [16J). It is however worthwhile recalling the 
proof: every stationary solution q of (jl .9f) satisfies 



l -q'{9) sm(9 - 9')q{9') dO'^j q(9) = C , (1.15) 

for some constant C . Since we know that q(-) > 0, then (|1.15p yields 

^ (log $(*))' -Jr(jf cob(0- 9')q(9')d9'J = (1.16) 

which implies C = 0. At this point, by playing on the rotation symmetry, we may assume 
that J s q{9) sin(0) d6 = 0, so that any stationary non-negative solution q(-) with prescribed 
first Fourier cosine coefficient f s q{9) cos(9) d9 equal to r satisfies 

-q(0)' - Kr cos{6)q{9) = 0. (1.17) 

A solution to (|1.17p is proportional to exp(2Kr cos(9). By normalizing (q(-) is a probability 
density) and recalling the constraint on the first Fourier cosine coefficient we get to (|1.13|) - 
(fTTIjh 

Remark 1.6. Remarkably, a generalization of (|1.13p - (|1.14p holds also in the disordered 
case [19J. The key to such a derivation, like in the step above, is in the identification 
of the order parameter r, that captures the degree of coherence (or synchronization) 
of the oscillators. In statistical mechanics terms this is nothing but the fact that the 
Hamiltonian Hjv(ip) may be rewritten as (1/2N) j Si ■ Sj = (l/2A r )(J^ i Si) 2 , with 
Si = (?R.exp(iipi),Q exp(iifi)). However, if one considers an /i-model (cf. Remark II. ip . 
the Hamiltonian cannot be expressed any longer as a function of the total magnetization 
Y^i Si- For the identification of the order parameter in this more general context we refer 
to [6J. 

1.3. The gradient flow viewpoint. For our purposes the following fact is of crucial 
importance: (jl.9p can be rewritten in the gradient form 



d t q t (9) = V 



(1.18) 



where we use V for dg for visual impact, SQ(q)/5q(9) is the standard L 2 Frechet derivative 
of the functional Q and 

HQ) ■= \ f q(9)logq(9)d9-^ [ cos(0 - 9')q{9)q{9') d9 d9' . (1.19) 
^ Js 2 Js 2 

Note that J- : L 2 (§>) — > R is Frechet differentiable at q(-) for example if q(-) is continuous 
and q(-) > and Proposition 13.11 guarantees that the evolution may be cast in the form 
(|1.18p for t > 0. A direct consequence of f 1 1 . 1 8 f) is that 



at 

A first consequence of this observation is: 
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Proposition 1.7. // there exists t2 > ii > such that q^i-) = 9t 2 ('); then there exists a 
constant C and a value r satisfying (I1.14|) such that qt(-) = q{- + c) for every t>t\. 

Proof. By Proposition 13.11 (guaranteeing smoothness and positivity of the solution) and 
by (|1.20p we see that ^(qt) is constant for t G [ii,t2], so that V5J r (qt)/5(qt(0)) = for 
every 9 and t G [ti,^]. But this implies that 

(log q t {9)) sin(0 - 9')q t {9') d9'^j = , (1.21) 

which is precisely (|1.16p with C = 0. Therefore for every t G [ii, £2] there exists a constant 
7(i) such that g t (0) = g(0 + 7(i)), with </(•) as in (fTTTH|> - (|TTTll ). Since q{- + j(tt)) is a 
stationary solution, the claim follows. □ 

Two observations are in order: 

(1) Proposition 11.71 generalizes to the non-disordered /i-model, when the latter is re- 
versible (see Section [4]), in the sense that the hypotheses imply V {5T / 5qt) = and 
this condition identifies all the stationary solutions. 

(2) Proposition 11.71 shows that there is no non-trivial stationary solution to (jl.ip when 
C; — £> £ a non-zero constant. This is simply because, otherwise, we would have 
a solution to fll.2|) of the form q(- — t£), with g(-) non-constant, which violates 
Proposition II .71 This is of interest also because it is not clear that Proposition [TTT] 
generalizes to disordered models. Clarifying the link between non-reversibility and 
coexistence of stationary and rotating solutions appears also to be an intriguing 
question. 

1.4. On synchronization stability. The main result that we present addresses the im- 
portant issue of the stability of the non-trivial stationary profiles q(-), more precisely of 
the stability of the invariant manifold {q(- + #o)}6» es- I n the literature we find a full 
analysis of incoherence stability [22] (also in presence of disorder) as well as an analysis of 
synchronized profiles as bifurcation from the incoherent 1/2-7T profile (we refer to pQ and 
the several references therein). Our aim is to have a detailed non-perturbative analysis of 
the linearized evolution operator in the non disordered case, for every K > K c = 1. 

To address such an issue we observe that the linearized evolution ut(-) around q(-) 
obeys the equation dtUt{9) = LxUt(9) with Lx a linear operator with domain -D(L^) := 
{C 2 (S; M) : / s u{9) &9 = 0} defined as 



q{9) f sin (9 - 9') u (9') d9' + u{9) [ sin (9 - 9') q (9') d0' 



Lq U (9) = -Au(9) + KV 

(1-22) 

It is easy to verify that Lqq' = 0, and this corresponds to the rotation invariance of the 
problem. However, what we are going to prove is that the remaining part of the spectrum 
is also real and it lies on the negative semi-axis. In order to make precise statements 
about Lq we introduce the Hilbert space j/g of distributions u such that u = W, 
with U G L 2 (S;M). Of course the derivative is taken in the sense of distributions and U 
is determined, given u, only up to a constant: we remove this uncertainty by stipulating 
that f s (U(9)/q(9)) d9 = 0. The norm of u G H_ 1A/q is defined by 

IM| 2 -i,i/, := f^d9, (1.23) 
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and the scalar product of u and v £ H-i i/g is going to be denoted by ((u,v)): it is of 
course equal to j s (U(8)V(9)/q(9)) d6, with definition of V in analogy with U. We will 
come back in the next section with more on H_ii/g, but what one can verify directly is 
that D(L^) and L^D(L^) are subsets of H_iij q and that Lg- is symmetric as an operator 
on H_i i/q, that is 

((u,L^v)) = ((v,L 4 u)), (1.24) 

for every u,v £ D(Lg) (for an explicit expression see (|2.14jl ). We will actually prove 
(Proposition I2.6| ) that L q is essentially self-adjoint. Moreover: 

Theorem 1.8. The spectrum of (the self-adjoint extension of) Lg is pure point and it lies 
in (— oo,0]. The value is in the spectrum, with one- dimensional eigenspace (spanned, as 
we have seen, by q' ) and the distance between zero and the rest of the spectrum is of at 
least 

X(K) = (l-^(l-O) (1 - (WKr)n 

K 2Kr 2 exp(8Kr) + exp(4Kr) (1 - (I Q (2Kr))- 2 ) ' K ' 1 

We stress that Theorem 11.81 holds as soon as there is a non-trivial solution r to (|1.14j) , 
that is for every K > 1. We have: 

X(K) K y^ and X(K) K ^ eM -^ + 2) . (1.26) 

Numerically increases till K = 1.033 . . ., where it reaches the value 0.0028 . . ., and then it 
decreases. 

The paper is organized as follows: in Section [2] we study Lg and prove a spectral gap 
inequality, the essential self-adjontness of the operator and the fact that the spectrum is 
pure point. The nonlinear evolution properties mentioned in this introduction are treated 
in Section [3] and the (ir)reversibility issues are considered in Section HI 

2. Synchronization stability 

In this section we prove the main result (Theorem ll.8j) . We assume K > 1 and, for 
simplicity, we drop the hat from q(-), so that a stationary solution is denoted by q(-). 

2.1. Some properties of the stationary profile. We first rewrite (jl.l3)l - (|1.14j) by using 
the Bessel function notation: 

#):=^^yexp(2^cos W ), (2.1) 

and r G (0, 1) is the unique positive solution of 

r := V(2Kr), with V(x) := (2.2) 

We have used the standard notation for the modified Bessel functions of order and 1, 
explicitly 

I v (x) := -!- / (cos(6»))^exp(xcos(6»)) d0, for v = 0, 1. (2.3) 

27r Jo 

As already mentioned before, uniqueness of r is a non-trivial fact that follows from [16, 
Lemma 4], that establishes in particular the concavity of $(•) on the positive semi-axis. 
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One can therefore define, via (|2.2|) ). the function [l,oo) 3 K ^ r(K) E [0,1) (one sets 
r(l) := by continuity). We have that, for K > 1, 1 — K(l— r 2 ) E (0, 1/2) or (equivalently) 

f^l < r(K) < f^T. (2.4) 

These bounds are easily checked for K close to 1 and K large, and and the numerical plots 
of the three functions appearing in (|2,4p ca be found in Figure [TJ We could not find quick 
proofs of (|2.4p : a proof of the upper bound is a byproduct of one of the arguments that 
we develop below (see the proof of Lemma l2.2p . while we prove here the lower bound by 
using Bessel functions properties. 



r{K) 




Figure 1. The function K i— > r(K) and (dashed lines) the bounds of (|2.4[l . The upper 
bound is proven, the lower bound is verified for K close to 1 and K large. 

Proof of (|2,4p . lower bound. By the change of variables (r, K) (r, 2Kr) =: (r, y) we see 
that what we have to prove is equivalent to showing that 



y>0 



m z {y) + -^{y)-l > 0, 

y 



Apply now the identity [23] 



(2.5) 



h{y) _y 

Uy) 2 I 1 + vhMl 



so that 



1 - ~¥(y) 

y 



1 _l 2 h(y) 
1 + yh{y) 



>0), 



and therefore (12. 5p is equivalent to 



V yhyy)J h{y)h 



(y) 



> i, 



(2.6) 



(2.7) 



(2.8) 



KURAMOTO MODEL AND MEAN FIELD SPIN XY DYNAMICS 



9 



where the intermediate step follows from the identity yl^iy) + 2ii(y) — yIo(y) = 
But this is equivalent to h(y) / Io(y) > I2 {y)/1 1 (y) for y > 0, a fact that is proven in [TH 
(3.8)]. □ 

In the sequel we will also use the notations 

J{9) := -Ksm(8) and J{6) := jfcos(0) , (2.9) 

so that 

Y = J*q, (2.10) 
and with this change of notation (|1.22p reads 

L q u = -u" -(q(J*u)+u(J* q))' . (2.11) 

2.2. Rigged Hilbert spaces and We now introduce a rigged Hilbert spaces 

structure [4, pp. 81-82]. The pivot (Hilbert) space is H := {u £ L 2 (S) : Ju = 0} (of 
course the scalar product is (u,v) := J uv and the norm is denoted by || • H2). The second 
Hilbert space we consider is V := H\ jq , closure of the set of periodic C 1 functions u such 
that J u = with respect to the norm 



IMIl, g := J j^ u '^ ( 2 - 12 ) 

so that V C H and the canonical injection of V into H is continuous (by the Poincare 
inequality). Note that V is dense in H. We consider then the dual space V' of V and the 
duality functional in V defined by (p u {v) := (u, v) for every given u € H (ip u : V — > M). 
We can define T : H — > V' by setting Tu := <p u . One can then show that T{H) is 
dense in V and T injects H into V in a continuous way [H p. 82]. This injection allows 
considering H as a subset of V, by identifying ti and Tu. Moreover if u G H we have that 
= ||Tu||y can be made more explicit: given u £ H we call IA the primitive of H 
such that J U/q = and we observe that 

\\u\\ v > = sup^^ = sup f = \ — 1 (2-13) 

vev \\v\\v veH hq Jfq( v >y \ J q 

where the last step follows on one hand by the Cauchy-Schwarz inequality (establishing 
the upper bound) and by choosing v' = U/q in the supremum (establishing the lower 
bound). 

As already mentioned in the introduction, the scalar product in V' = H—11/q is denoted 
by ((•,•)). We observe also that these steps allow the precise identification of the functions 
in H—i Uq' u € H_i ij q if and only if u = U' (in the sense of distributions), with U E L 2 (S). 

At this point it is crucial to observe that D(L q ) (recall that D(L q ) is the subset of 
periodic C 2 functions u such that j u = 0) is dense in H_ 11 / q and that for u, v S D(L q ) C 
H-11/q, we have (use (|2.10p ) 

((v, L q u)) = ((L q v,u)) =-lj™ + (v,J*u) . (2.14) 



In words: L q is a symmetric operator on H_i y 
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2.3. Estimates on the Dirichlet form. It is going to be useful to introduce also the 
Hilbert space L 2 ,^, which coincides with L 2 (S) as a set of functions, but we equip it with 
the scalar product 

(u,u) : = (u,u/q). (2.15) 

We therefore introduce, for u G D(L q ), the Dirichlet form V(u) := —((u,L q u)). By (12.140 
we have 

V{u) = ^ (u, u) - (u, J* u) . (2.16) 

Our aim is to bound from below T>(u) and we start with two technical lemmas. The 
first one yields the spectral decomposition of (u, J * u), viewed as a quadratic form on 
L 2 . 

1/9 



Lemma 2.1. We have the orthogonal decomposition 

2 

where 



L\/ q = Vo e x V 1/2 ® x Vkr-i/2, (2-17) 



V := <8^a + ^2(aj cos(j6) + bj sin(j0)) : ^ o| + b) < oo > , (2.18) 
[ i>2 j ) 

and both V\/2 and Vr-i/2 are one dimensional subspaces generated respectively by 8 ^ 
s'm(8)q(8) (= —q'(8)/2Kr) and by 8 i— > cos(8)q(8). Moreover, when u £ V\ we have 

J*u = -u. (2.19) 

Proof. The L^-orthogonality statements Po_LT4/2 an d Vo-LVk-i/2 follow directly from the 
orthogonality in L 2 of the family {cos(j6>), sin(j^)}j = o i i v ... Instead ^i/2-L^k-i/2 because 
JZ 7T q(8)cos(8)sm(8)d8 = 0. 

The validity of (|2.19p follows by direct computation: for u € Vo 

J * u{8) = K cos{8) cos(8')u(8')dd' + Ksm{8) / sm{8'))u(8') d8' , (2.20) 
Jo Jo 

which is equal to zero because u does not contain the first harmonics. The other two cases 
follow by using the same trigonometric identity and the following two (clearly equivalent) 
identities: 

J 2 J q(8) sm 2 (8) d8 = ± ^ q(8) cos 2 (0) dd = (l - ^\ . (2.21) 

□ 

The second lemma is more technical and its interest will become clear in the proof of 
Proposition 12.31 

Lemma 2.2. We have 

2K — 1 

min (1 + u, 1 + u) = (1 + u, 1 + u) = (2tt) 2 — - ^— , (2.22) 
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where 

2tt ( 1 —) ?Trr 

fi W = - 1 " f 2 71 i \ \ 9(0) + q(9) 008(0), (2.23) 

Proof. We have to minimize a quadratic functional under the linear constraints (l,tt) = 
and ii E Vo. This corresponds to the three constraints: 

/ i 2-7r /•2ir /*27r 

i u(9)d6 = 0, / cos(9)u(9)d6 = and / sm(9)u(9)d9 = 0. (2.24) 
J ■/ o Vo 

The extrema (minima, by convexity) of such a problem can be found by the Lagrange 
multipliers method and they are of the form 

u{Q) = -1 + Xq(9) + fiq(6) cos(0) + Tjq(d) sin(0), (2.25) 

with A, (j, and r\ three real numbers. The constraints (|2.24p . via (|2.2p and (|2.2ip . yield 
the linear system r\ = 0, A + [ir = 2ir and 2KXr + [i(2K — 1) = 0, which has a solution 
if and only if 2ii~(l — r 2 ) — 1 7^ 0. Since the minimum exists for every K and since 
2K (1 - r 2 ) - 1 -> 1 for K \ 1 we see that 2.ff (1 - r 2 ) - 1 > for every K (this is the 
upper bound in (|2.4|) ). The proof is completed by making A and \i explicit and using that 
the expression in (12321) is equal to A 2 + 2A/xr + ,u 2 (l - (1/2*0). □ 

The following is one of our main statements: 

Proposition 2.3. There exists ck £ (0,1/2) such that, if u £ L\/ q is such that (l,u) = 
J Q 27r u = 0, then 

V(u) > c K (u - u 1/2 ,u - Ui/ 2 ) , (2.26) 
where n 1 / 2 := q' (u,q'} / {q',q')- In particular, V(u) > 0. 

Proof. Lemma [2. II shows that, if we write u = (according to the decomposition 

(|2.17p : of course Vij 2 = ^1/2) we have 

V(u) = -(K-l)(v,v) + ±(v ,v ). (2.27) 

We write v(9) = cq(9) cos(#) and vq = a + ^2j >2 (cij cos (i^) + &j srn (j#)) and, with this 
notations, one directly sees, using the definitions (|2.ip and (|2.2p . that the constraint 
(1, u) = is equivalent to 

re = -2vra . (2.28) 

Let us observe that we can assume c / 0: if c = then v — V\i 2 = vq and, in view of 
(|2.27p . (|2.26p holds. From now on we perform estimates for arbitrary, but fixed, values of 
c 7^ (hence is fixed too). 

Note then that (v,v) = — (1/(2K))), by (|2.21|) . and we can therefore write 

1 \ ~2 , °n / u 2 , , u 2 



W>- + *<l + %l + %), (2.29) 

where of course we have set u 2 '■= vq — a®. By Lemma 12.21 we therefore obtain that 

min V(u) = V(v + a Q (l + u)) , (2.30) 

u:(l,u)=0 
given c 
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where we recall that c and ao are related via (|2.28|) . Note that, by (|2.4|) . ci(K) > 
for K > 1. For sake of compactness let us introduce vo := ao(l + Vo), so that 

u = (vo — vo) + vo + 1^/2 + v and, looking back at (|2.27|) . we see that 

V(u) = -(K -l)(v,v) + -(v ,v ) + (vo,vo-v ) + - (v -vq,vq-v ) 

= V{v + vo) + (vo, v - vo) + - («o -vq,v q - vo) , / 2 .3l) 

= min V{u) + (u 0) «o - v Q ) + - (v - v , v - %) , 

u:(l,u)=0 2 
given c 

which implies in particular 

(v , v - v ) + ~ («o - «o, «o - wo) > , (2.32) 



Since a lengthy computation yields 
V(y + v ) 



l-K{l-r 2 ) =:c K £ (0,1/2), (2.33) 



(v,v) + (ucwo) 
where ck S (0, 1/2) is just a restatement of (|2.4|) . we get 

V{u) = C K ((v,v) + (v ,Vo)) +(v ,V -Vo) + - {vq - V ,V - V ) 

> ck ((v,v) + (vo,vq)) + 2ck (vo,v - Vq) + ck (vq - v ,v - vo) (2.34) 

= C K ((V,V) + (VQ,V )) , 

and the proof is complete. □ 

For the next result we point out that, by the definition of in terms of rigged Hilbert 

spaces, we know that there exists c > such that ((it, it)) < c(it, it) for every u E I? with 
(1, it) = 0. Therefore 

((«,«)> < Cp (u,u) , (2.35) 
with cp = cmaxg. This can be proven directly and cp can be made explicit, in fact 
J U/q = tells us that U(9 ) = for some 6> , so that ((it, it)) < (mine/) -1 J e ° +27{ W 1 (we 
are looking at IA as a periodic function with domain R) and the Poincare inequality tells 
us that Jq° +27T M 2 is smaller than J^ +2lT ipi') 2 = (u, it). Therefore we can choose 

cp = = exp UKr). 2.36 

mint? 



Lemma 2.4. For u G £ 2 (§) sitc/t t/iai (1, it) = ice /laue 

1/9 ((?',?')) 

where C > is given by 



(2.37) 

#-1,1/9 
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Of course Proposition 12.31 and Lemma 12.41 yields the spectral gap inequality: for every 
u such that ((it, </')) = we have 

V(u) > c K C 2 ((u,u)). (2.39) 

Proof. Set e = q' / (q' , q') 1 , so that u%/ 2 = (u, e) e. By (|2.35p we see that (|2.37p follows if 
one can show 



u M 

((e,e)) e 



1/9 



(2.40) 



ll n ~ u 1/2\\l 2 , > C P C 

1 1/9 

= c P C ]j ( u - u i/2>u ~ Wi /2 > + (V e) ^ ^ ^ , 
and this is equivalent (note that cpC 6 (0, 1)) to 

< U - u 1/2 , u - « 1/2 > > Co ((„, «) - I;!) 2 _ C„ ( «"-^- e)> ) 2 , (2.41) 

with Co := {cpC) 2 /{\ — (cpC) 2 ). By the Cauchy-Schwarz inequality and by (|2.35|) we 
have 

((it - n 1/2 ,e)) 2 < ((u-u 1/2 ,u-u 1/2 )}((e,e)} < tip (u - u 1/2 , u - u x/2 ) ((e, e)), (2.42) 

so that (I2.41|) holds if C c P /((e,e)) < 1, which is equivalent to CpC 2 < ((e,e))/(cp + 
((e, e))). This is a condition on C that we can verify explicitly by using ((e,e)) = (1 — 
I 2 (2Kr))/{2Kr 2 ). □ 

2.4. Self-adjointness and spectral properties of L q . Let us recall that H is the space 
L 2 with zero average constraint and that L q is viewed as an operator on H_ii/ q D ff. 
The first result is a technical lemma: 

Lemma 2.5. Fix K > 1. There exists c > suc/t i/iarf /or every u £ H , we have 

((«,«)) > cK /2 || 2 , (2.43) 
where U\/ 2 is the orthogonal projection, in L> 2 ^ q , on V\j 2 . 

Proof. Using the explicit expression for ui/ 2 (Proposition I2.3p . and 

j{\ogq)'u = -J (log q)"U, (2.44) 



in, a') 



we see that 

IK/2II2 < (maxg) (u 1/2 ,uy 2 ) 



^ f f (log q)"u) < 2vrmaxg 2 max |(log g )-f f =; 1 



U 2 /q =: -((«,«)), (2.45) 



where the last step is the definition of c. □ 
Proposition 2.6. L g zs essentially self-adjoint. 
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Proof. We start by introducing for u, v G D{L q ) 

St(v,u) := ((v,(l-L g )u)) = - j M + l -j V ^--jvJ*u. (2.46) 

The right-most expression shows that £\(u, v) is well defined as long as u,v G H (this 
generalizes the definition of £i(-, •) and we will take this definition from now on). Moreover 
£i(-, •) is a continuous and coercive bilinear form on H x H, that is there exists c G (0, 1) 
such that 

£i(u,v) < -IMhIMb; and £i(u,u) > c||ii|| § . (2.47) 

The second inequality follows from Proposition 12.31 ad Lemma 12.51 Now observe that for 
every / G H_ii/ q the linear form v \— > ({v, /)), from H to M, is continuous (it is continuous 
also as a map from H—i t i/q to M) and therefore, by the Lax-Milgram Theorem [4j Cor. V.8], 
we have that there exists a unique u G H such that 

£ x (v,u) = ((v,f)}, for every v G H. (2.48) 

Since we can write 

((v,f)) = - J dfl, (2.49) 
from (I2T46D . (I2TI91) and (|2^8|) we see that 



« 2 9 («) 

for (Lebesgue) almost every 9. Since u £ H, the primitive oiU/q is C 1 and its (weak) 
second derivative is square integrable. If / G H(c H_i i/ q ) then the same is true for the 

right-hand side in (|2.50p . Since J* u is C°°, we see that u is C 1 (more precisely, has a C 1 
version). So the left-most term in (I2.50P is at least C 3 . If now we assume that / is C°, we 
can therefore conclude that u G C 2 . 

To sum up: if / is periodic, C , with J / = 0, then n G D(L q ) and (1 — L g )u = /, which 
follows by taking applying dg(q(8)dg ■ ) to both terms in (|2.50p (and by using (I2.10p ). Since 
such functions / are dense in H_i ij q , we see that the range of 1 — L q is dense, so that its 
kernel is {0}, and this implies that L q is essentially self-adjoint ([U Prop. VII.6]). □ 

Proposition 2.7. The spectrum of L q is pure point. 

Proof. We are going to prove this by showing that the resolvent of L q is compact, namely 
that (A — Lq^ 1 is compact for A in the resolvent set. It suffices to prove such a result for 
one value of A |12l p. 187] and we choose A = 1, which is in the resolvent set thanks to 
Proposition 12.31 and Proposition 12.61 So let us consider / := (1 — L g ) _1 u, u G H_i i/ q , so 
that / is in the domain of 1 — L q and we have 

((/, (1 - L q )f)) = «/,u». (2.51) 

But, by (I2T4TD . ((/, (1 - L q )f)) is bounded below by c||/|||, so that 

cUh - urn- c-mt s c*«' {Zb2) 

where we have used the continuous injection of H into V' = H_i i/ q (C is the constant 
arising when comparing the norms of these two spaces). Therefore (1 — Lq) -1 maps 
sequences that are bounded in H_ 11 / q to sequences that are bounded in H. We are 
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therefore left with showing that the embedding of H into H_i i/q ls compact. This just 
follows by the Cauchy-Schwarz inequality: for every v G H we have \V(0) — V{9')\ < 
IMI2 Vr — ®'\ an( ^ ' since we are on a bounded interval with periodic boundary conditions 
and J V/q = 0, this yields that {V : v G H and \\v\\2 < const.} is a compact subset of C° 
(Ascoli-Arzela Theorem), and hence of Ly q . That is, a bounded subset of H is a relatively 
compact subset of H_i i/ q and this completes the proof. □ 

3. The nonlinear evolution 
We need the following result on the nonlinear evolution: 

Proposition 3.1. For every uq G .A/fi(§) there is a unique element v. o/C°([0, 00); A4i (§)) 
sitc/i i/tat (|1.8j) is satisfied for every F G C 2 (§) and every t > 0. Moreover for t > t/ie 
measure vt is absolutely continuous with respect to the Lebesgue measure and, if we denote 
by Qt(') its density, the function (t,6) 1— > qt{9), with domain of definition (0, 00) x S, is 
C°°, strictly positive and it solves ()1.9|) . 



Proof. The uniqueness result is proven for example in [15J (with § replaced by R) and 
in [8] (for a general, bounded or unbounded, domain subset of M. d ). The case of periodic 
boundary conditions is not treated explicitly but the argument of proof goes through 
with minor modifications. Existence is established by the tightness result on the particle 
system, which follows by methods that are by now standard: it suffices in fact to show 
that, for every (smoooth) F, {J s F dvN,t}N is tight (see e.g. [T3l [T5l E] ) and this follows 
immediately from the fact that the drift in (jl.2p is bounded. 

For the regularity and positivity aspects we use the fact that, if for t > there is a 
solution to 

dtu(t,6) = l J^l + K ^- e [{H{t,e))u(t,9)\, (3.1) 

with 

H(t,0) := / sin(0-0> t (d0') = s in(0) f cos(9')v t (d9 f )-cos(9) [ sm(9')v t (d9'), (3.2) 
(note that H(-,-) is continuous in time and C°° in space) such that for every F £ C°(S) 

lim / F(9)u(t,9)d9 = [ Fdu , (3.3) 
*\° J s J 

then one directly verifies that v. G C°([0, T\;Mi(§)), defined by v t (d6) = u(t,9)d9 for 
t > and vq = is a solution to (jl.8p . Hence v. = v. by uniqueness. 

Equation (13. If) is a parabolic linear partial differential equation on which there is much 
literature. For our purpose the results by D. G. Aronson in [2] turn out to be particularly 
relevant. In particular, Aronson shows that the equation (13.1j) (consider 9 G M. for the 
moment) admits a fundamental solution T(t, 9; s, 9') (t > s) so that the solution at time t 
can be written as j R T(t, 9; 0, 9')u(0, 9') d9', at least when u(0, •) G Lf oc (M) (of course the 
solution has to be interpreted in a weak sense: see [21 pp. 608-609] for the very general set 
up in which such a result it is proven). A number of results on the fundamental solution 
are proven [21 Section 7] and notably that it is a continuous function in (t, 9) for t > s and 
that, if we assume that sup te r s s+T ]^ g]R if(t, 9) =: M < 00, there exists C = C(T,M) > 
such that 

' e^(- C ^ 2 ) < TM;^) < ° a JJ±Zm , (3.4) 



Cy/t^S V Vt^S J y/t^S V C^ft 
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for every 9, 6' and every t G (s,s + T]. Moreover ([21 Corollary 12.1]) the result applies 
also to initial data that are measures: namely, if p is a measure on R such that for t > 
the map 9 i— ► J" K F(t, 9; 0, 9')p,( d9) is a function in Lf oc (M) then such a map defines a weak 
solution to (|3,ip on (0,T] xR (namely a weak solution to (|3.ip on [to,? 1 ] x R for every 
*o€(0,T)). 

Let us now specialize to our case: H(t, •) is smooth and 27r-periodic, so that T(t, 9; s, 9') = 
T(t, 9 + 2ir; s, 9' + 2ft). Let us apply the result we have just stated with p defined by re- 
quiring that its restriction to [2irj,2ir(j + 1)) coincides with the image of the measure vq 
under the application [0, 2ir) 3 9^9 + 2nj G [2ft j, 2ft(j + 1)), for every j. Therefore for 
t > 

/ T(t,9;0,x)p(dx) = V / T(t,9;0,9' + 2ftj)u(d9') =: v(t,9), (3.5) 



and, by (|3.4p . v(t, •) is bounded as soon as t > 0. Therefore «(-, •) is a weak solution and, 
in turns, v(t,9) d9 coincides with vt(d9), which therefore has a representation in terms of 
the fundamental solution T and this implies not only that the solution becomes a bounded 
function as soon as t > 0, but also (by the lower bound in A3.4B ) that it is strictly positive. 

At this point, since we know that the solution is bounded, the smoothness in both 
variables of the solution (for t > 0) may be derived by standard methods: this issue is 
taken up for example in [3] for a slightly different evolution equation, or, more generally, 
in P2j. □ 

4. On the irreversibility of the Kuramoto model 

In this section we consider the Kuramoto fa-model, i.e. (Il.ip with general drift, as in 
Remark 1 and a = 1. 

Existence of a unique invariant probability measure (for each fixed N) is a well known 
fact, but one can actually prove that such an invariant measure has a positive C°° density 
p : S N — ► (0, oo). These issues are treated in detail for example in [14, Ch.s 3 and 
5], where one finds also an extensive treatment of the entropy productions for Markov 
processes (with references to the vast literature on the subject). The entropy production 
rate e p for a stationary process X is defined as the limit as T — > oo, when it exists, of 
the relative entropy of the law of {Xt} t ^\o,T] with respect to the law of {XT-t}te[o,T} > 
divided by T. For a large class of models e p takes the form of the steady state average 
of time integral of the square of a suitable flux. This is true also in our case, namely 

e p = \ J&n Ejli J j (f) 2 P(<£) where 

2 N d 
Jj (p) := 2ij -j^Yl h ^3 ~ <Pi) ~ qT; lo g P (^) • 

i=i J 

The key point is that e p = if and only if the system is reversible \14:\ Th. 5.4.6] (of course 
reversibility calls for specifying an invariant probability with respect to which the system 
is reversible, but in our set-up there is only one invariant measure). Therefore our system 
is reversible if and only if Jj(-) = for every j: let us spell it out 

d 2 N 

— log/? (ip) = 2£j Kfj ~ <Pi) for ever y 3 and <£• ( 4 - 2 ) 

^ i=i 

This expression directly implies that J^h(9)d9 = 2ir£j for every j, that is £j does not 
depend on j. Without loss of generality we may therefore assume = for every j (recall 
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Remark 1 1.2p . which entails j s h(9) d6 = and therefore the primitive h of h is 2-7r-periodic 
(we make the arbitrary choice h(0) = 0). By integrating (|4.2p we obtain that 

2 „ 

log P (ip) = j^Yl h (w ~ + c j fe) ' ( 4 - 3 ) 

j 

where Cj (93) does not depend on tpj. For j = 1 and ifi fixed for i = 3, 4, ... we can rewrite 
as 

2 ~ 

logp(vi,^2) = JtHvi ~ V2) +51 (yi) +52(^2) , (4.4) 

where 51(^1) := (2/iV) Yli>3 h&j ~ fi) an d 92(^2) '■= c\ (</>)■ We can of course repeat the 
same steps with j = 2 obtaining thus 

2 ~ 

log ^1,^2) = ^(^2 - yi) + /i(yi) + /2(<^2) , (4.5) 

with /1 and /2 defined in analogy with g\, g 2 (but we are simply interested in the fact 
that they are smooth functions from § to R). From (14. 4|) and (14. 5h we infer that 

/i((^i - ip 2 ) - h(ip 2 - ipi) = f(tpi) + g(<p2), (4.6) 

for suitable smooth functions / and g from S to R. This tells us in particular that 
f(c + 6) + g(c) does not depend on c (it is equal to h(8) — h{—9)) and, therefore, that 
f{x) + g(x) is a constant and f(c + x) — f(x) = /(c) — /(0) for every c and x (that is / is 
constant, since it is continuous and periodic). We have therefore reached the conclusion 
that i — > h(8) — h{—9) is a constant, which has therefore to be zero. 

We sum up the argument we have just developed in the following statement: 

Proposition 4.1. For every N the dynamics defined by (ll.l|) . generalized as in Re- 
mark [Ll\ is reversible if and only if the following two conditions are satisfied: 

(1) fl = £7 for every j; 

(2) h(-) — £1 : R ^ R is an odd function. 
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